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@) TECHNIQUE POUR L'ESTIMATION DE PLUIE A PARTIR D'UN RADAR METEOROLOGIQUE A DIVERSITE DE 
POLARISATION. 

(§) Procede pour {'estimation d'un taux precipitant au 

moyen d'un radar caracte rise en ce que, connaissant la pa- et on determine R2R connaissant K2R. 
rametre de proportionnalite entre I'attenuation specifique K 
= A2D et le taux de variation differential de la phase, on de- 
termine une estimation de la valeur K2R20 & partir de 
I'equation 
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on determine une estimation de K2R a partir de I'equa- 
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La presente invention est relative aux techniques meteorologiques 
pour ('estimation d'un taux precipitant au moyen d'un radar. 

Pour des raisons de cout, on envisage aujourd'hui d'utiliser des 
radars a longueurs d'ondes courtes (bande X) permettant de reduire la taille 
5 d'antenne pour une resolution angulaire donnee. Malheureusement, a ces 
frequences, les ondes radar sont attenuees par la pluie ce qui rend la 
restitution du taux precipitant plus delicate qu'en bande C (ou S) ou les 
effets d'attenuation sont beaucoup plus limites. 

Le probleme de la correction d'attenuation a ete envisage par 
10 Hitschfled et Bordan qui en ont donne des 1954 une solution 
mathematiquement rigoureuse, sous Phypothese d'une relation 
parametrique entre la reflectivite Z et ('attenuation specifique K, du type loi 
puissance : K = a Zo b . 

On pourra a cet egard avantageusement se referer a leur publication 
15 ; HITSCHFELD - BORDAN, 1954, "Errors inherent in the radar 
measurement of rainfall at attenuating wavelengths", J. Meteor, 11, pp. 58- 
67. 

Malheureusement, dans une application reelle ou les donnees sont 
bruitees et ou Ton peut etre en presence d'une erreur d'etalonnage (rnerne 
20 faible), cette solution se revele numeriquement instable. 

Pour pallier cet inconvenient, il a deja ete propose de fonctionner en 
diversite de polarisation et d'exploiter la reflectivite differentielle ZDR, qui est 
le rapport des reflectivites en polarisation H et V. 

Au cours de sa chute, une goutte de pluie subit un aplatissement 
25 aerodynamique qui a pour consequence que sa section efficace de diffusion 
est plus grande en polarisation H qu'en polarisation V. Dans un radar 
operant a frequence attenuee, le parametre ZDR resulte de deux effets de 
tendances opposees : d'une part, I'effet de retrodiffusion differentielle 
(tendant a accroitre ZDR avec la taille moyenne des gouttes), et d'autre part 
30 I'effet d'attenuation differentielle (de tendance opposee). 

Pour un exemple recent de determination d'un taux precipitant a 
partir d'un radar a double polarisation, on pourra avantageusement se 
r eferer a la demande de brevet FR-2 742 876, ainsi qu'a la publication : 
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"Polarimetric Radar at Attenuated Wavelenghts as a Hydrological sensor" - 
M. SAUVAGEOT - Journal of etmosyheric and oceanic technology - vol. 1 3 
-p. 630-637, 1996. 

Toutefois, la mesure du coefficient ZDR n'est pas aisee. 
5 La dynamique de variation de ZDR est d'environ 1 a 2 et couvre done 

quelques dB tout au plus. En outre, cette mesure requiert un inter- 
etalonnage des voies H et V au dizieme de dB ; Integration sur de 
nombreux echantillons independants pour temperer la fluctuation statistique 
du signal, et de tres bonnes performances de I'antenne (lobes secondaires 
10 a 30 dB - aller - en dessous du lobe principal). 

En alternative, il a egalement ete propose, notamment par la 
publication ; 

"Differential propagation phase shift and rainfall rate estimation" - M. 
SACHINANDA, D.S. ZRNIC - Radio science - vol. 21, n°2, p. 235-247, 

15 March-April 1986. 

d'exploiter un autre parametre, a savoir la phase differentielle (notee <Ddp). 
En effet, I'aplatissement de la goutte de pluie n'affecte pas seulement la 
section efficace de retrodiffusion induisant la dissymetrie o H > ov (ou a H et 
ov represented les sections efficaces en polar H et V, respectivement), 

20 mais aussi la propagation de I'onde radar dont le vecteur d'onde est affecte 
tant pour sa partie imaginaire (I'attenuation specifique est plus forte en 
polarisation H qu'en polarisation V, effet d'attenuation differentielle 
precedemment mentionne), que pour sa partie reelle (ou la dissymetrie 
entre H et V se traduit par une variation differentielle de la phase le long du 

25 trajet. Or, la derivee de Odp le long du rayon (taux de variation de Odp. 
notee Kdp) est quasi-proportionnelle au taux precipitant (et done Odp, au 
taux precipitant integre), ce qui explique Tinteret de le mesurer. 

L'estimateur de Kdp presente de nombreux avantages : il est 
insensible a rattenuation le long du trajet ; il est insensible a une erreur 

30 d'etalonnage du radar ; il est beaucoup moins affecte que ZDR par les lobes 
secondaires de I'antenne ; il n'est aucunement affecte par un masque partiel 
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du faisceau d'antenne (ce qui arrive systematiquement quand on opere a 
faible elevation). 

II presente cependant ('inconvenient majeur d'etre fortement bruite 
car il resulte de la derivation le long du faisceau de cpdp. La mesure de Kdp 
5 requiert done un long temps d'integration, incompatible avec I'application 
hydrologique qui impose un balayage rapide du faisceau radar. 

L'invention propose quant a elle une technique qui permet de pallier 
les inconvenients des techniques anterieures et qui met en oeuvre un 
traitement simple, fiable et robuste, en particulier vis-a-vis du bruit 
10 statistique de mesure (ce qui permet un balayage rapide). 

Plus particulierement, l'invention propose un procede pour 
('estimation d'un taux precipitant au moyen d'un radar bipolaire, caracterise 
par les differentes etapes suivantes : 

- on mesure au moyen dudit radar bipolaire, sur un intervalle [ri, r 0 ] donne 
15 de rayon r de trajet par rapport audit radar, la phase differentielle Q>dp et la 

reflectivite attenuee Z selon au moins I'une des polarisations H ou V, 

- on determine une estimation de la valeur K(r 0 ) de ('attenuation specifique 
en r 0 a partir de Tequation 

avec 7(5, r) = 0A6b£ Z b (u)du 

ou Z(r) est la reflectivite mesuree le long de la radiale et b est I'exposant de 
25 la loi puissance [K(r 0 ) = a Z(r 0 ) b ] liant I'attenuation specifique en r 0 a la 
reflectivite non attenuee Z(ro), 

et ou y est le parametre de proportionnalite entre I'attenuation specifique et 
le taux de variation differential de la phase, 

- on determine une estimation de K(r) a partir de I'equation 



30 
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£j2-*3&>., fr .,.> 

Kir) K{r,) 

- on determine le taux precipitant R (r) connaissant K(r). 
Avantageusement, on mesure la reflectivite attenuee Z pour Tune et 

5 I'autre des deux polarisations H ou V et on determine pour Tune et I'autre de 
ces deux polarisations le taux precipitant R(r) pour un rayon de trajet donne 

r. 

Uinvention concerne egalement un dispositif pour I'estimation d'un 
taux precipitant comportant un radar bipolaire, ainsi que des moyens de 

10 traitement, caracterise en ce que ledit radar comporte des moyens pour la 
mesure de la phase differentielle Odp et de la reflectivite attenuee Z selon 
au moins Tune des polarisations H ou V et en ce que les moyens de 
traitement mettent en ceuvre les differentes etapes de traitement du 
procede selon la revendication 1 . 

15 D'autres caracteristiques et avantages de invention ressortiront encore 
de la description qui suit, laquelle est purement illustrative et non limitative 
et doit etre lue en regard des desstns annexes sur lesquels : 
- la figure 1 illustre un dispositif pour la mise en oeuvre du procede conforme 
a I'invention ; 

20 - les figures 2a a 2c, 3a a 3c, 4a a 4c sont des graphes, dans lesquels on a 
suppose que le faisceau radar intercepte deux cellules de pluie de forme 
gaussienne sur lesquels on a porte les parametres Z, K et R pour 
differentes valeurs du parametre N 0 (defini plus loin dans la description) de 
distribution dimensionnelle des gouttes de pluie. 

25 - la figure 5 est un graphe sur lequel on a porte une simulation de 

variations pour le coefficient N 0 en fonction du rayon de trajet r 

- les figures 6a a 6c sont des graphes sur lesquels on a porte, pour le 
cas ou N 0 est variable, les valeurs obtenues pour les parametres Z, K et R ; 

- les figures 7a a 7c et les figures 8a a 8c sont des graphes sur 
30 lesquels on a porte les parametres Z, K et R obtenus pour des biSis absolus 
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de 2dB et 5 dB respectivement et pour une valeur de N 0 et une prise egale 
a 8000 ; 

- les figures 9a a 9c et 10a a 10c sont des figures semblables a 
celles des figures 7a a 7c et 8a a 8c dans le cas de biais relatifs pris 

5 respectivement egaux a 0.2 dB et 0.5 dB ; 

- les figures 1 1a a 1 1c et 12a a 12c sont des graphes sur lesquels on 
a porte les valeurs obtenues pour les parametres Z, K et R dans le cas de 
restitution a parttr d'un nombre d'echantillons egal a 60 et d'un nombre 
d'echantillons egal a 10 ; 

10 - les figures 13a a 13c sont des graphes sur lesquels on a porte les 

valeurs obtenues pour les parametres Z, K et R dans le cas d'un taux de 
pluie maximum de 200 mm/heure et d'un signal bruite. 

Le procede qui va maintenant etre decrit met en ceuvre un traitement 
sur la phase differentielle Odp et la reflectivite attenuee Z H (ou Z v ) selon la 
15 polarisation H (ou V) (reflectivite Z dans la suite du texte). 

Dans une premiere etape, on determine pour la distance radiale r a 
partir de I'antenne radar un intervalle [r 1t r 0 ], dans lequel les donnees Z et 
cDdp sont disponibles continuellement. 

On sait que ('equation d'Hitschfeld - Bordan s'ecrit : 

20 

u" , ^* = _o.46aZ* m 
dr W 

avec u = Z/Zo ou Zo est la reflectivite attenuee et Z la reflectivite apparente, 
et ou a et b sont les deux parametres de la loi de puissance K = a Zo b , dans 
25 laquelle K est Pattenuation specifique. On rappelle que a est fortement 
variable, tandis que b est stable. 

La solution de cette equation (1) est, lorsque Ton integre entre deux 
bornes n et r 2 : 

=-bA6ab^Z b cls (2) 

30 

En exprimant (2) entre r et r 0( il vient : 
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K(r) K(r 0 ) 

ou 

5 /(r, r 0 ) = 0.466 J'°Z 6 ds (4) 

et ou K = aZo b 

^equation (3) fournit un estimateur de K(r) que Ton va chercher a 
contraindre par reformation de <Ddp en ajustant K(r 0 ). Pour cela on 
transforme (3) en 

10 

Z b (r) 

K(r) = K(r 0 ) ^2 (5) 

Or, la derivee de la phase Odp par rapport au trajet r, ou taux de 
15 variation de la phase differentielle Kdp, est une fonction quasi-lineaire de 
I'attenuatioft specifique K (c'est-a-dire K H ou K v , selon la polarisation sur 
laquelle oh travaille). 

Par consequent, 

20 J^^)A = y(^(^o)-M r .)) < 6 > 

ou y est un parametre de proportionnalite connu 
soit : 

Jr *Z(r 0 ) + K(r 0 )J(s,r Q ) 

25 

A partir de cette estimation de K(r 0 ), on determine K(r) en tout point 
de I'intervalle [n, r 0 ] a partir de I'equation (3). On remarque que Tinteret de 
cette formulation est d'une part, qu'elle fournit une solution K(r) 
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independante du parametre a qui est precisement le plus fluctuant, et 
d'autre part, qu'elle fait appel a une mesure differentielle de la phase 
differentielle, beaucoup plus aisee a atteindre qu'une rnesure absolue. 

Ceci permet d'estimer la valeur du taux precipitant R(r), en utilisant 
5 par exemple a cet effet des abaques donnant R(r) en fonction de K(r) pour 
plusieurs valeurs du parametre a, lequel peut etre lui-meme determine a 
partir de la relation : 

K(r 0 ) = aZ(r 0 ) b 

En mettant en ceuvre le traitement qui vient d'etre decrit pour Tune et 
10 I'autre des deux polarisations H et V, on dispose, in fine, de deux 
estimations de R qu'il est possible de croiser pour valider une valeur finale. 

Le traitement qui vient d'etre decrit a ete teste avec un appareillage 
du type de celui illustre sur la figure 1 qui comprend : 

- un emetteur 1 ; 

15 - un diviseur de puissance 2 a 3dB recevant le signal en sortie de 

I'emetteur 1 ; 

- un extracteur de mode (ou orthomode) 3 relie par deux voies H et 
V aux sorties d'un diviseur de puissance 2 ; 

- une antenne 4 a cornet transmettant ou recevant les deux 
20 polarisations H et V ; 

- des circulateurs 5 (te de filtrage) interposes, sur les voies H et V, 
entre I'extracteur de mode 3 et le diviseur de puissance 2 ; 

- un dephaseur 6 interpose, sur la voie V, entre le diviseur de 
puissance 2 et I'extracteur de mode 3 ; 

25 - des moyens de reception 7 relies, pour les voies H et V, aux 

circulateurs 5 ; 

- un oscillateur 9 qui synchronise I'emetteur 1 et les moyens de 
reception 9 ; 

- des moyens 8 d'echantillonnage, puis de traitement des sorties I 
30 et Q de ces moyens de reception 7, ces moyens de traitement 8 

permettant la determination a partir des parametres ZH, ZV et 
4>DP, lesquels verifient en notation complexe: 
Z h< v = [Iii,v 2 + Qh,v 2 ] 1/2 et ODP = 1 arg[(l H + iQH).(lv-iQv)] 
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ou l H , v et Q Hl v sont les valeurs obtenues par echantillonnage des 
sorties I et Q des voies H et V. 

Ce sont ces moyens de traitement 8 qui mettent en oeuvre le 
traitement de determination qui a ete decrit ci-dessus. 
5 Les parametres de fonctionnement de ce dispositif etaient les 

suivants : 

- seui! de detection du radar a 0 dB SNR et a 1 0 km : 0 dBZ ; 

- nombre d'echantillons independants Nj dans la cellule de 
resolution : 1 0 ou 60 

10 - resolution radiale : 75 metres 

- bruit sur la mesure de <t> 0 p : ± 1 degre (pour Nj=60), et ± 3 degres 
(pour Ni=10) 

- distribution dimensionnelle des gouttes exponentielle : N(D) = N 0 
exp (-AD) [ou D est le diametre moyen des gouttes et ou diametre 

15 A est un parametre], avec N 0 constant le long de la radiale (avec 

les valeurs suivantes : 800, 8000 [valeur de Marshall et Palmer], 
80000, ou continument variable proportionnellement a R. 
Dans ces conditions, les differents resultats suivants ont ete 
obtenus : 
20 1) N o constant . 

Les figures 2a a 2c, 3a a 3c, 4a a 4c presentent les taux de pluie 
restitues pour un No constant sur la radiale et valant respectivement 800, 
8000, 80000. Dans ces 3 cas, on voit que les taux de pluie restitues sont 
tres proches des taux de pluie simules, et que la valeur du parametre N 0 
25 restituee est a chaque fois relativement proche de celle utilisee pour les 
simulations. 

2) Nn variant avec le taux de pluie. 

Les inventeurs ont egalement teste le traitement sur des simulations 
30 toujours realistes mais realisees avec un No dependant du taux de pluie, ce 
qui est assez conforme a la realite mais qui rend M assez perilleuse H la 
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restitution du taux precipitant. La relation utilisee entre N 0 et le taux 
precipitant est la suivante : 

5 1000 

Pour des taux de pluie variant de 0 a 80 mm/h, on voit que le 
parametre No explore une tres large gamme, couvrant plus d'une decade. A 
titre dedication, on a represents sur la figure 5 les variations de N 0 sur la 
10 radiale simulee. Les restitutions obtenues avec le traitement precedemment 
decrit sont presentees sur les figures 6a a 6c. L'ecart maximum observe est 
de I'ordre de 20%, ce qui semble raisonnable dans le cas de simulations 
"aussi defavorables" 

15 3) Effets des biais relatif et absolu. 

Les figures 7a a 7c et 8 a 8c presentent les resultats obtenus avec 
un biais absolu de 2dB et de 5dB respectivement. Une valeur de biais 
absolu de 2dB induit une erreur non negligeable sur I'estimation du 
parametre N 0l mais ne penalise pas la restitution, puisque les taux de pluie 

20 restitues ne different de ceux qui ont ete simules que de quelques 
pourcents. En revanche, un biais absolu de 5dB semble redhibitoire. Les 
resultats obtenus dans le cas de biais relatifs de 0.2 dB et de 0.5 dB et qui 
sont presentes sur les figures 9a a 9c et 10 a 10c, montrent la faible 
sensibilite de cet algorithme a une erreur d'etalonnage relative entre les 

25 deux polarisations 

4) Effets du bruit . 

Les inventeurs ont teste I'effet du bruit sur les restitutions dans 
le cas d'un nombre d'echantillons de 60 (figures 11a a 11c) et dans le cas 
30 plus defavorable de 10 echantillons (figures 12a a 12c). Dans ces deux cas, 
la mesure la mesure de la phase integree presente un bruit de 3 degres. La 
restitution est realisee tant que le rapport signal a bruit est suffisant. Pour un 
taux de pluie maximum de 80 mm/h, c'est le cas jusqu'a une distance 
radiale d'environ 19 km. Les resultats obtenus confirment le comportement 
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tres stable de ralgorithme KDP. Dans ces conditions des coefficients de 
correlation superieurs a 0.95 et des biais inferieurs a 4 mm/h sont tout a fait 
satisfaisants. 

5 5) Cas de tres forts taux de pluie . 

Les inventeurs ont verifie la robustesse de ralgorithme dans le cas 
de forts taux de pluie, atteignant 200 mm/h et un signal bruite. Dans ce cas, 
le signal se noie assez rapidement dans le bruit, et evidernment il devient 
impossible de restituer la deuxieme cellule de precipitations a cause de la 
10 tres forte attenuation du signal radar par la premiere. Cependant les taux de 
pluie obtenus jusqu'a environ 10 km sont coherents et Tensemble de la 
premiere cellule precipitante est tres convenablement restitue (figures 13a a 
13c). 

Le traitement a egalement ete teste dans d'autres bandes que la 
15 bande X et notamment en bande C. Les resultats obtenus sont tout aussi 
satisfaisants qu'en bande X. 

Ainsi, avec le traitement qui vient d'etre decrit, on dispose d'un 
estimateur du taux precipitant qui est particulierement stable vis-a-vis du 
bruit et/ou du biais sur les reflectivites et sur la phase. Ce traitement 
20 parvient, y compris dans les cas difficiles (fort taux de pluie, N 0 variant 
rapidement avec le taux de pluie) a restituer le taux de pluie de maniere tres 
encourageante. 

Le traitement des donnees est possible aussi longtemps que le 
rapport signal a bruit est superieur a 1. Dans ces conditions, le biais sur 

25 I'estimateur de R n'excede pas 5%, et son ecart type 10%. La performance 
de I'estimateur est independante de N 0 (quand on le suppose constant le 
long de la radiale), et n'est que legerement deterioree dans les conditions 
de N 0 variable. Un etalonnage absolu du radar a ± 2 ou 3 dBZ suffit pour 
assurer ces performances. Cet etalonnage n'est requis qui pour affiner la 

30 relation R-K a partir de laquelle on procede a I'estimee de R. Aucune 
specification n'est requise en ce qui conceme I'etalonnage relatif des voies 
HetV. 
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L'estimateur de taux de pluie qui vient d'etre decrit est 
particulierement avantageux dans le cas de restimation de la pluie sur les 
bassins versants pour la surveillance des crues et la gestion des ressources 
en eau. Cette surveillance et cette gestion reposent aujourd'hui sur des 
modeles de ruissellement/infiltration dont la principale limitation quant a la 
performance predictive tient 3 une mauvaise connaissance de la lame d'eau 
entrant dans le systeme. 
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REVINDICATIONS 

1. Procede pour ('estimation d'un taux precipitant au moyen d'un 
radar bipolaire, caracterise par les differentes etapes suivantes : 

- on mesure au moyen dudit radar bipolaire, sur un intervalle [r 1( r 0 ] donne 
de rayon r de trajet par rapport audit radar, la phase differentiate d>dp et la 
reflectivite attenuee Z selon au moins Tune des polarisations H ou V, 

- on determine une estimation de la valeur K(r 0 ) de Tattenuation en r 0 a 
partir de ('equation 

ou: I(s t r 0 ) = 0A6bj'°Z b (u)du et b est Texposant de la relation 

puissance (r 0 ) en r 0 

K(r 0 )=aZ(r 0 ) b 
ou K(r) est Tattenuation specifique en r 0 . 

et ou y est le parametre de proportionnalite entre Tattenuation specifique et 
le taux de variation differential de la phase, 

- on determine une estimation K(r) de Tattenuation specifique en fonction de 
r a partir de Tequation 

K{r) K(r 0 ) oJ 

- on determine le taux precipitant R (r) connaissant K(r). 

2. Procede selon la revendication 1 , caracterise en ce qu'on mesure 
la reflectivite attenuee Z pour Tune et Tautre des deux polarisations H ou V 
et en ce qu'on determine pour Tune et Tautre de ces deux polarisations le 
taux precipitant R(r) pour un rayon de trajet donne r. c 

3. Dispositif pour ('estimation d'un taux precipitant comportant un 
radar bipolaire, ainsi que des moyens de traitement, caracterise en ce que 
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ledit radar comporte des moyens pour la mesure de la phase differentielle 
0>dp et de la reflectivite attenuee Z selon au moins Tune des polarisations H 
ou V et en ce que les moyens de traitement mettent en oeuvre les 
differentes etapes de traitement du precede selon la revendication 1 . 
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